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AI-7Si-Mg 合金 凝固 过 程 形 核 模 型 建立 及 
校 晶 生长 过 程 数值 模拟 
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摘 要 针对 铝 合金 砂 型 铸造 较 低 冷 速 特点 , 通过 实测 和 分 析 不 同 凝固 条 件 下 的 冷却 曲线 , 建立 了 适用 于 铝 合金 形 核 密度 随 
最 大 形 核 过 冷 度 呈 指数 性 变化 的 形 核 函数 . 通过 与 Pandat 软 件 热力 学 动力 学 `. 平 衡 相 图 数据 库 相 耦合 , 并 利用 空间 坐标 变 
化 等 算法 , 建立 了 适用 于 三 元 铝 合金 二 维 、 三 维 枝 晶 生长 的 CA 模型 . 在 该 模型 中 , 同时 考虑 了 溶质 扩散 、 成 分 过 冷 、 曲 率 过 
冷 . 唱 体 择 优 取向 以 及 不 同 组 元 之 间 相 互 作 用 等 重要 因素 的 影响 . 利用 建立 的 形 核 和 生长 模型 , 模拟 了 Al-7Si-0.36Mg 合金 
在 不 同 凝固 条 件 下 的 二 维 枝 晶 演化 及 形 貌 特征 , 描述 了 溶质 组 元 的 分 布 特征 以 及 定量 地 预测 了 二 次 枝 晶 臂 间 距 的 变化 , 并 
与 实验 结果 进行 了 对 比 . 三 维 枝 晶 的 模拟 结果 有 效 反映 了 枝 晶 空间 结构 复杂 性 和 多 样 性 , 并 与 实验 结果 吻合 良好 . 

关键 词 三 元 铝 合金 , 形 核 模型 , 元 胞 自动 机 , 枝 晶 生 长 , 二 次 枝 晶 臂 间距 
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ABSTRACT Dove to the extensive applications in automotive and aerospace industries of Al-7Si-Mg casting al- 
loys, its understanding of the dendrite microstructural formation is of great importance to control the desirable mi- 
crostructure and thereby to modify the performance of castings. In this work, through analyzing the measured cool- 
ing curves in different cooling conditions of A1-781-0.36Mg ternary alloy during sand casting, a theoretical nucle- 
ation model correlated maximum nucleation undercooling with the nucleation density is proposed. Besides, a 2D 
and 3D cellular automaton (CA) model allowing for the quantitatively predicting dendrite growth of ternary alloys 
is presented. This model introduces a new tracking neighboring rule algorithm to eliminate the effect of mesh de- 
pendency on dendrite growth. The thermodynamic and kinetic data needed in the simulations is obtained by cou- 
pling with Pandat software package in combination with thermodynamic/kinetic/equilibrium phase diagram calcula- 
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tion databases. This model has also taken account the multi-component diffusion, constitutional undercooling, cur- 
vature undercooling, dendrite preferential growth angles as well as the effect of interactions between the alloying 
elements etc. This model is applied to quantitatively simulate the dendrite growth with various crystallographic ori- 
entations of Al-78i-0.36Mg ternary alloy in 2D and 3D during polycrystalline solidification, and the predicted sec- 
ondary dendrite arm spacing (SDAS) shows a reasonable agreement with the experimental results. The experimen- 
tal observed complicated and diverse dendrite microstructure occurring in solidification process can be well repro- 
duced by this 3D-CA model which has considered the effects of various preferred growth orientations, the interac- 
tions of adjacent dendrites as well as the influence of S/L interface anisotropies. The simulated results effectively 
demonstrated the abilities of this model in prediction of dendrite microstructure in ternary alloys. 

KEY WORDS ternary aluminum alloy, nucleation model, cellular automaton, dendrite growth, secondary den- 


drite arm spacing 


金 因 其 轻 质 、 高 强 、 耐 腐蚀 、 易 成 型 等 一 系 热 扩 散 、 溶 质 扩 散 以 及 校 晶 生 长 动力 学 等 , 建立 较 


列 优点 , 目前 作为 一 种 功能 结构 材料 , 广泛 应 用 在 ”复杂 的 描述 枝 品 固 / 液 界面 形态 以 及 溶质 分 布 的 数 


航空 航天 汽车 家电、 食品 等 行业 上 . Al-7Si-Mg (i 学 方程 , 从 而 在 更 大 尺度 上 实现 对 枝 唱 形 貌 演化 过 
量 分 数 , %, 下 同 ) 合 金 ( 如 ZL114, A356, A357) 属 于 程 的 模拟 , 更 加 接近 于 工程 应 用 实际 . 因此 , 近 些 年 
典型 的 ALSi 系 铸造 铝 合金 , 因为 其 含有 较 高 的 Si ”来 发 展 较 快 , 在 实际 中 得 到 了 广泛 的 应 用 李斌 


ru 


, 使 液态 金属 具有 较 好 的 流动 性 , 可 以 成 型 大 等 "在 大 量 实验 统计 分 析 基 础 之 上 , E RE 


型 薄 壁 复杂 的 铸件 , 同时 由 于 微量 Mg (0.3%~0.7%) — XEZLIIAA 铝 合金 瞬时 形 核 模型 和 相应 的 CA 法 , 并 
元 素 的 加 入 , 可 通过 形成 MgsSi 析 出 相 来 进行 时 效 。 利用 阶梯 件 开 展 了 模拟 和 实验 研究 , 两 者 吻合 良 
强化 处 理 , 是 汽车 缸 体 、 动 车 齿轮 箱 等 重要 零 部 件 好 . 霍 亮 等 扩 和 吴 备 武 等 思 分 别针 对 镁 合金 hcp 特 


的 典型 材料 . 近年 来 , 随 着 汽车 工业 的 不 断 发 展 ,对 点, 通过 定义 不 同 的 邻居 单元 类 型 和 捕获 规则 , 建 
Al-7Si-Mg 合 金 铸件 的 性 能 提出 了 更 高 的 要 求 , 因 — 立 了 适用 于 镁 合金 hcp 结 构 的 枝 品 生长 CA 模型 , 并 
此 , 关于 该 合金 在 不 同 工 艺 下 的 凝固 行为 、 微 观 组 A 建立 的 模型 分 别 模拟 了 不 同 冷却 速率 下 铂 合 
织 演化 以 及 性 能 等 方面 一 直 备 受 关注 金 的 凝固 组 织 以 及 不 同 压 铸 工 艺 下 镁 合金 的 枢 唱 
铸件 力学 性 能 主要 受 凝 固 组 织 的 影响 , 而 凝固 。 E, 并 与 实验 结果 进行 对 比 . 石 玉 峰 等 "以 透明 合 
组 织 又 决定 于 合金 成 分 及 铸造 工艺 参数 . 为 了 提高 。 金 NHCI-H,O 水 溶液 定向 凝固 实验 结果 为 基础 , 通 


合金 的 强化 效果 , 实际 合金 一 般 都 含有 多 种 合金 元 。 过 建立 一 种 改进 的 CA 模型 模拟 和 分 析 了 不 同 温 度 


素 , 所 


以 研究 不 同 凝固 条 件 对 多 元 合金 凝固 组 织 的 ”梯度 、 校 唱 择 优 取向 以 及 不 同 抽 拉 速度 对 柱状 唱 生 


影响 规律 具有 重要 的 理论 意义 和 工程 价值 . 对 于 长 形态 的 影响 . 然而 , 实际 应 用 的 合金 多 数 是 三 元 


Al-7Si-Mg 合金 来 说 , 该 合金 的 初生 相 凝 固 温 度 


z 


或 是 多 元 合金 , 而 上 述 的 模型 都 主要 局 限 在 二 元 合 


间 相 对 较 罕 , 在 共 品 形成 之 前 , 熔 体 中 的 校 晶 组 织 金枝 唱 的 模拟 , 所 以 建立 多 组 元 合金 校 晶 生长 模 
没有 足够 的 时 间 进 行 枝 晶 臂 的 充分 生长 , 所 以 该 合 W, 不 管 是 从 理论 分 析 还 是 工程 应 用 等 角度 都 是 非 
金 的 校 晶 辟 组 织 一 般 比 较 发 达 . 枝 品 的 一 二 次 术 唱 。” 常 必要 的 . Zhu 等 中 通过 把 CA 模型 与 热力 学 相 平衡 


SN 


E, 尤其 是 二 次 梳 唱 臂 间距 , 是 用 于 表征 凝固 计算 数据 库 (PanEngine) 相 耦合 , 建立 了 适用 于 描述 


昌 织 特征 的 重要 参数 , 其 大 小 直接 影响 着 成 分 .第 ”多 元 铝 合金 校 唱 生长 和 微观 偏 析 的 数学 模型 , 并 与 


2 从 而 影响 铸件 的 使 用 性 能 Yr ”Scheil 模 型 预测 结果 进行 对 比 . 但 是 该 模型 引入 的 动 
年 来 , 二 次 校 蝇 辟 间距 与 铸件 的 力学 性 能 之 间 的 天 ”力学 生长 系数 是 经 验 性 参数 , 该 值 选 取 的 大 小 会 影 
系 引 起 了 人 们 的 广泛 关注 响 计 算 结 果 . 文献 [16~18] 在 计算 多 元 合金 校 唱 生长 


随 着 计算 机 技术 的 快速 发 展 , 数值 模拟 技术 作 W, 主要 基于 界面 局 部 平衡 求解 界面 生长 速度 , 从 而 


快速 发 展 . 通过 对 


1 一 种 预测 凝固 组 织 演变 的 有 效 途 径 , 得 到 了 ”来 计算 固 相 分 数 的 变化 . 但 该 方法 需要 求解 界面 法 
不 同 凝固 条 件 下 的 微观 组 织 演变 H, Wie N A E 特别 是 对 三 


和 最 终 组 织 特征 进行 定量 模拟 , 从 而 对 铸件 局 部 点 “，” 维 杭 唱 的 计算 , 同时 该 方法 一 般 需 要 通过 从 代 的 方 
力学 性 能 进行 预测 已 经 成 为 可 能 . 目前 , 微观 组 织 。 式 进行 求解 . 因此 , 现 有 的 多 元 合 4 HE CA M 


数 


LM 


直 模拟 的 主要 方法 有 元 胞 自动 机 (CA) 法 和 相 场 ” 型 还 需 进一步 的 改进 和 完善 . 实际 的 凝固 过 程 涉及 


(PF) 法 , 其 中 , CA 法 上 共有 较 高 的 计算 效率 , 通过 耦合 ”到 三 维 复杂 的 多 枝 品 空间 结构 , 所 以 建立 能 够 描述 
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三 维 多 枝 晶 生长 过 程 的 CA 模型 , 使 得 计算 结果 更 加 
符合 实际 凝固 组 织 , 是 目前 枝 晶 模拟 的 一 大 难点 . 此 
外 , 针对 实际 合金 体系 , 结合 工业 生产 实际 开展 模拟 
研究 , 增强 模型 的 实用 性 也 是 今后 的 发 展 方向 号 . 

形 核 是 凝固 过 程 的 第 一 步 , 决定 了 相 的 析出 次 
序 和 具体 分 布 , 是 准确 预测 凝固 组 织 特征 的 先决 条 
件 , 因此 建立 合适 的 形 核 模 型 显得 尤为 重要 . 以 往 
凝固 过 程 微观 组 织 模拟 大 多 采用 基于 Gauss 分布 的 
连续 形 核 模型 3, Oldfield"" 提 出 的 用 于 模拟 激 冷 灰 
铸铁 的 连续 形 核 模 型 或 是 Stefanescu 等 外 提出 的 瞬 
时 形 核 模型 . Gauss 连续 形 核 模型 的 晶 粒 密度 n(A7) 
可 以 表示 为 过 冷 度 AT 的 函数 , 但 在 利用 该 模型 确定 
唱 粒 密度 或 是 最 终 的 形 核 数 目 时 , 首先 必须 测定 平 
均 形 核 过 冷 度 AT,,, 过 冷 度 标准 偏差 AT, UJ o 
形 核 质点 密度 ns, 而 这 3 个 参数 比较 难以 测 得 , 需 
要 测定 不 同 过 冷 度 AT 下 的 形 核 密度 进行 拟 合 得 到 ， 
而 且 形 核 密度 计算 结果 受 这 3 个 参数 的 偏差 比较 敏 
感 , 以 至 于 在 微小 的 枝 唱 计算 区 域 (mm? RUE) WW, 该 
模型 计算 得 到 的 最 终 形 核 数目 往往 与 实际 偏差 较 
大 . Oldfield 形 核 模型 是 根据 实验 数据 统计 获得 , 预 
测 结果 相对 较为 准确 , 但 是 该 模型 中 的 形 核 参 数 缺 
乏 明 确 的 物理 意义 . 
本 工作 首先 基于 实验 数据 的 统计 , 建立 适用 于 
Al-7Si-0.36Mg 合金 在 较 低 冷却 速率 条 件 下 ( 砂 型 铸 
造 条 件 ) 凝 固 的 形 核 模 型 , 并 与 CA 枝 唱 生长 模型 相 
耦合 , 使 得 模拟 过 程 更 加 符合 凝固 过 程 物理 本 质 . 
通过 与 Pandat 软件 热力 学 动力 学 数据 库 相 耦合 ， 
建立 适用 于 模拟 三 元 铝 合金 梳 唱 生长 的 二 维 、 三 维 
CA 模型 , 并 通过 空间 坐标 变换 等 算法 , 以 消除 网 格 
各 向 异性 对 枝 晶 生长 的 影响 , 实现 任意 角度 枝 晶 生 
长 过 程 的 描述 . 利用 建立 的 形 核 和 生长 模型 模拟 


模型 中 的 参数 缺乏 明确 物理 意义 . 本 工作 借鉴 该 模 
型 思想 , 建立 了 适用 于 铝 合金 等 轴 唱 凝固 的 形 核 模 
型 , 同时 模型 中 的 参数 具有 明确 的 物理 意义 , 并 且 
能 够 方便 地 耦合 到 枝 唱 生长 模型 当中 . 
为 了 建立 经 过 Al-Ti-B 细 化 处 理 的 铝 合金 熔 体 
凝固 过 程 的 形 核 模型 , 简化 了 复杂 的 形 核 本 质 , 提 
出 如 下 的 假设 条 件 : (1) 不 考虑 晶 核 的 溶解 与 合并 ， 
即 假设 原始 形 核 数 目 等 于 最 终 的 唱 粒 数目 , 同时 忽 
略 细 化 剂 的 衰退 影响 ; (2) 根据 文献 [23,24], 经 过 Al- 
Ti-B 孕育 的 铝 合 金 熔 体 , ga-Al 的 初始 形 核 过 冷 度 非 
常 小 , 基本 小 于 0.5 K, 甚至 可 以 达到 0.01 K, 因此 可 
以 认为 当 温 度 一 旦 低 于 液 相 线 温度 时 , 形 核 就 开 
始 , 而 当 温 度 达 到 7 时 , 形 核 结束 ; (3) 实验 中 观察 
到 TiB; 的 形 貌 呈现 六 边 形 板 状 , 假设 熔 体 中 TiB; 是 
唯一 形 核 核心 “9, 其 {0001} 基 面 如 图 1a 所 示 . 
许多 学 者 “3 针对 经 过 Al-Ti-B 孕育 的 铝 合 金 熔 
体 中 的 ga-Al 形 核 过 程 进行 了 相关 的 实验 研究 , 研究 
表明 , ga-Al 只 在 TiB, 的 {0001} 基 面 上 依附 形 核 . 首 
先 , 形 核 的 核心 会 沿 着 TiB; 衬 底 向 周围 生长 , 直至 
履 盖 整个 {0001} 基 面 , 之 后 的 生长 主要 通过 减 小 国 
液 界 面 的 曲率 半径 ns 来 实现 , 即 图 1b 中 的 高 度 不 
断 增 大 . 随 着 h 增 大 , rs 减 小 ,但 是 根据 自由 生长 理 
论 中 ,在 固定 的 过 冷 度 AT 下 , rs 不 能 小 于 临界 核心 
半径 7 BUramr.rRmEATX 
2ys. 
T NSQAT Q) 
式 中 ，ys 表示 固 / 液 界面 能 ，AS、 fe HR DPI 
4t. Ni. 通过 式 (2) 可 以 看 到 , 只 有 通过 增 大 AT 减 小 
r, Q-Al 才 能 得 到 进一步 生长 . 当 w-Al 在 TiB, 衬 底 上 
达到 如 图 lb 中 所 示 的 半球 形状 时 , rs 达到 最 小 值 ， 
即 ra —70.5d. (H P, 四 为 图 lb 所 示 的 形 核 衬 底 直 径 )， 


T 


Al-7Si-0.36Mg 合金 在 不 同 凝固 条 件 下 的 凝固 组 织 ， 
并 与 实验 结果 进行 比较 . 
1 形 核 模型 建立 

Oldfield2 基 于 实验 统计 , 建立 了 预测 铸铁 和 球 
铁 凝 固 过 程 中 共 晶 团 或 是 石墨 球形 核 密度 的 函数 ， 
该 模型 认为 , 形 核 密 度 W 可 表示 为 最 大 形 核 过 冷 度 
AT, 的 函数 ; 


N=y(AT,) (1) 
式 中 , y 和 是 通过 实验 统计 获得 的 形 核 常数 ; A7,= 
TuTon 其 中 , Ti 为 平衡 液 相 线 温度 , Tu 对 应 于 准 
却 曲 线 上 的 最 低 形 核 温度 . 由 于 该 模型 基于 实验 数 
据 拟 合 获 得 , 计算 简单 , 结果 相对 较为 准确 , 工程 应 
用 性 强 , 目前 广泛 应 用 在 各 种 模拟 软件 当中 , 但 该 


之 后 晶体 就 能 自由 生长 . 根据 假设 (2), 在 形 核 阶段 
当 温度 降低 至 Fu 时 , 形 核 过 冷 度 达到 最 大 值 AT 
此 时 能 够 形 核 的 衬 底 直 径 最 小 , 可 表示 为 : 

d= 6) 


从 式 (3) 可 知 , 车 已 知 TiB; 形 核 衬 底 的 尺寸 分 布 , 就 
可 以 建立 形 核 密度 和 AT, 之 间 的 函数 关系 . 
Greer 等 中 基于 扫 面 电镜 (SEM) 图 像 分 析 技 术 ， 
对 Al-5Ti-1B (质量 分 数 , %) 孕 育 剂 中 的 TiB; 颗 粒 的 
尺寸 分 布 进行 了 大 量 的 数据 统计 , 并 用 指数 函数 
行 数据 拟 合 , 得 到 了 如 下 的 分 布 函数 : 
No step (4) 


式 中 , 四 表示 尺寸 分 布 的 特征 直径 , 为 单位 体积 炊 
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(a) (b) 


Minimum-radius 


hemispherical cap 


1 TiB, 有 效 形 核 面 及 ga-Al 依 附 在 TiB; 衬 底 上 的 形 核 过 程 示 意图 


Fig.1 Schematic of active circular nucleant surface of TiB; (a) and the growth of a- Al phase from such a nucleant site 


with a diameter of d (b) (L and S represent liquid and solid phases, respectively. / is the height of the nucleus on 


the nucleant site) 


体 中 总 的 TiB, 颗 粒 数 目 , Md) 73 El S ELTE TE dw~dx+ 
dd 之 间 的 TiB; 颗 粒 密度 . 函数 N(qy) 相 对 dd\ 的 曲线 
关系 如 图 2a 所 示 , 曲线 Nd\) 在 颗粒 直径 dcs. 
(其 中 , dws 和 ds 分 别 为 最 小 和 最 大 形 核 讨 底 直径 ) 
之 间 所 围 成 的 区 域 面积 表示 单位 熔 体 中 总 的 形 核 
衬 底 数目 Nu, 而 灰色 区 域 dy ~qdyws 之 间 所 围 成 的 区 域 
面积 则 表示 能 够 形 核 的 有 效 衬 底数 目 . 根据 假设 
(2), 当 熔 体 中 的 过 冷 度 达到 最 大 值 A 思 时 , 如 图 2b 
所 示 , 形 核 结束 , 此 时 颗粒 尺寸 dd\ 的 TiB; 颗 粒 
都 能 达到 形 核 条 件 , 而 尺寸 四 < d, 的 颗粒 达 不 到 所 
需 的 形 核 过 冷 度 . 因此 , 熔 体 中 只 有 少数 部 分 的 TiB， 
颗粒 能 够 作为 有 效 的 形 核 衬 底 . 铬 在 固定 的 孕育 条 
件 下 , 通过 降低 Ta 增 大 AT 则 熔 体 中 能 够 形 核 的 
有 效 衬 底 增多 , 合金 细 化 效果 增强 . 
通过 积分 , 可 以 得 到 有 效 的 形 核 衬 底数 目 , 从 而 可 了 AAA 
单位 熔 体 中 的 形 核 数 N 其 计算 公式 如 下 所 示 : Fig.2 Schematics of the size distribution of TiB, nucleus 


sites (a) and the corresponding cooling curve (b) 


Nd,) € 


N,min 
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N v7 In N (d )ddy (5) (dumin and dy, o —minimum and maximum nucleus 


EARBS), Ny 可 以 表示 为 : sites, d, —nucleus sites corresponding to the maxi- 
a mum nucleation undercooling, 7;,—liquidus temper- 
Ny =N, exp- ) (6) 


A T 


ature, 7,4—minimum nucleation temperature, 7— 


temperature, t—time 
4z Aya P ) 


- (7) 
AS. d - 
"NN ve 化 . 前 一 种 主要 受 形 核 条 件 的 影响 , 而 后 者 取决 于 
可 以 看 到 , 式 (中 的 2 个 参数 加 和 a 都 有 明确 。 生长 过 程 . 因此 本 工作 在 建立 形 核 模型 的 基础 上 
的 物理 意义 ,参数 4 由 ys AS URKI AMEN E GA KEMIA Aat BERE 
共同 决定 由 于 缺少 N AS 和民 参 数 的 数据 通过 采用 枝 品 生长 CA 模型 来 描述 枝 唱 在 不 同 凝 
E ' We M Md ， ”条 件 下 的 生长 过 程 并 预测 最 终 的 枝 晶 形 貌 . 本 模型 
无 法 获得 式 (6) 中 参数 的 具体 值 , 因此 本 工作 主要 利 主要 基干 以 下 合理 假设 来 对 枝 量 的 符 长 过 程 进行 
用 实验 数据 NN, 和 AT,, 通过 数据 拟 合 的 方式 来 求 得 MU CAN QU Ge 
aea 描述 : (1) 由 于 微观 计算 域 很 小 , 在 mm 尺度 , 因此 可 
T E 以 考虑 溶质 扩散 在 微观 计算 域内 满足 溶质 守恒 , 即 
A Aee E ENA 采用 封闭 边界 条 件 ; 同时 假设 整个 计算 域内 温度 志 
金 凝 固 组 织 细 化 程度 主要 体现 在 两 方面 ， 均匀, 而 温度 值 根据 实际 的 冷却 曲线 来 获得 ; (2) 由 
SAL, 另 一 方面 是 二 次 枝 晶 臂 间距 的 细 ”于 Al 和 Si 元 素 的 密度 差异 很 小 , 自然 对 流 对 枝 晶 生 
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长 影响 较 小 , 因此 不 考虑 液 相 流动 的 影响 ; (3) 忽略 
动力 学 过 冷 度 , 因此 模型 只 适用 于 模拟 低 Pkclet 数 
条 件 下 的 枝 晶 生长 . 
2. 多 组 元 溶质 扩散 方程 

微观 CA 单元 有 液态 、. 固 态 和 固 / 液 界面 3 种 状 
态 . 由 于 界面 单元 固 相 率 的 增加 排出 每 种 溶质 元 
A, 使 得 界面 前 沿 存在 溶质 梯度 , 从 而 推动 溶质 扩 
散 的 进行 , 因此 需要 对 每 种 元 素 分 别 进行 溶质 扩散 
方程 的 求解 , 其 公式 为 : 


ow? 


xxt 2 D; Vw?) wi - k) : (8) 


过 去 


AP, w PETE (si 或 是 Mg) 在 meo 
相 ) 当 中 的 溶质 含量 ; D; 73 p 4E P BS TP ORB E; 大 为 
固 相 分 数 ; ! 为 时 间 ; k 为 平衡 溶质 分 配 系数 , 多 元 
情况 下 它 不 再 是 常数 , 而 是 与 成 分 和 温度 有 关 , 需 
要 耦合 平衡 相 图 数据 库 来 获得 . 
2.2 生长 动力 学 及 固 相 分 数 增 量 的 计算 

对 于 界面 单元 , 成 分 过 冷 和 界面 曲率 的 共同 作 
用 可 以 使 界面 单元 满足 局 部 热力 学 平衡 条 件 , 因而 
界面 平衡 温度 可 以 表示 为 : 

T(i) - Tj ws. wt) - AT; - AT, (9) 
RP, TL wy) 为 平衡 液 相 线 温度 , 可 以 根据 Si 
和 Mg 的 成 分 值 I of , 从 平衡 相 图 上 读 取 ; TU 
是 计算 域 中 单元 的 实际 温度 值 , 通过 冷却 曲线 来 获 
得 ; AT, LAT, 分 别 对 应 于 成 分 过 冷 和 曲率 过 冷 ， 
它们 的 计算 公式 为 : 
AT, e= Ty us wg.) - Ty Qus wg.) ~ 

m Wy.) 


C 


X. — s (wr - w) (10) 
AT, = au y(n)+ Ae (11) 


RP, oT, (we wy )/ow; 这 一 项 表示 组 元 i 在 成 分 为 
w 时 的 液 相 线 斜率 ; wr 是 界面 单元 的 局 部 平衡 溶 
质 成 分 ; k 为 局 部 界面 平均 曲率 , 该 参数 的 计算 可 
以 参见 文献 [27,28]; 9 和 6 是 用 于 表示 界面 单元 空 


间 法 线 方向 的 角度 EE: / pod 


y oos yan) 


为 各 向 异性 函数 , 由 下 式 "给 出 : 
y(n) - AS. - 32): 


| 十 t (cos*0, + sin'0,(1 — 2sin'6, cos'0 D (12) 


1 
其 中 , TN Gibbs-Thomson 系数 ，s 为 界面 能 各 向 异 
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性 系数 . 
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界面 单元 回 相 分 数 的 计算 主要 有 两 种 方式 , 一 


种 是 基于 固 / 液 界面 推 


移 过 程 中 溶质 守 


面 生 长 速度 , 从 而 得 到 


谍 界 面 局 部 平衡 , 对 了 


局 部 杠杆 定律 


恒 方 程 求解 界 
固 相 率 的 增加 量 ; 另 一 种 是 考 
于 平衡 溶质 分 配 系数 小 于 1 的 合 
金 , 认为 界面 平衡 溶质 成 分 和 界面 实际 溶质 成 分 的 
差 值 是 枝 晶 生长 的 驱动 力 , 从 而 利 
直接 求 得 固 相 分 数 的 增 量 , 避免 了 求解 枝 晶 生长 速 
RE, 计算 效率 高 . 本 模型 主要 采 
界面 单元 在 某 一 时 间 步 长 内 由 于 固 相 分 
加 , 会 排出 部 分 的 Si 和 Mg 元 素 , 满足 如 下 方程 : 

ws -wg =w- ks)Afs 

wy, — Whg = y (E — ky )Afs 


j 后 一 种 方法 . 若 革 
数 的 增 


(13) 
(14) 


式 中 , A 为 界面 单元 在 每 一 时 间 步 长 的 固 相 分 数 增 


这 3 个 未 知 数 . 


量 . 联 立 式 (9)~(14), 可 以 唯一 求解 wt ， 


2.3 任意 角度 枝 晶 生长 描述 
唱 的 优先 生长 方向 是 随机 的 ， 


凝固 过 程 中 , 枝 


wy, 以 及 A 


如 图 3a 中 所 示 , 但 是 由 于 CA 网 格 的 各 向 异性 , 会 导 


致 枝 唱 的 生长 方 其 固有 的 优先 生长 方向 . 


了 消除 网 格 各 向 异性 


给 枝 品 生长 带 来 的 固有 影响 ， 


本 工作 采用 了 空 


= 间 坐 标 变换 和 人 
意 角度 梳 晶 生长 的 模拟 . 


i 心 算法 来 实现 任 


如 图 3a 所 示 , 对 于 
局 部 坐标 系 (xo Yo, Zo), 
而 世界 坐标 系 (x, yz) 沿 着 网 格 划分 方向 , 因此 2 个 
坐标 系 之 间 的 位 置 关系 可 以 通过 如 下 的 空 


人 


E 


为 


每 一 个 校 唱 都 赋予 其 各 自 的 


坐标 轴 治 着 各 自 的 [001] 方 向 ， 


旋转 矩阵 Rotate, Rotate», Rotate; 来 实现 . 


Rotate, = 


z 


其 中 , a, B, yz VA A 


所 示 . 在 二 维 条 件 下 ， 
可 以 简化 为 : 


o)" 


Rotate, = 


Rotate, = 


[8] A s 
cosy -siny 0 
m y cosy 0 
0 0 1 
cosB 0 sing 

| 0 1 0 (15) 
-sinB 0 cosB 


1 0 0 
0 cosa -singa 


0 sina cosa 


Zo 


KERA Z IRIR Euler £8, 如 图 


x Xo 
y |= Rotate," Rotate, Rotate;| yo (16) 


3b 


即 a=0, 8-0, li] A^ s 5 TEK R 


cosy -siny \ Xo 
cos y fs qu 
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Color online 


3 不 同 优先 生长 角度 枝 蝇 示意 图 及 局 部 坐标 系 (xo, yo, zo) 与 世界 坐标 系 (x, yz) 之 间 的 位 置 关系 


Fig.3 Schematic of the various growth orientations of different grains (a), and the location relationship between the local 


coordinate system (Xo, yo, zo) and the world coordinate system (x, y, z) (b) (a, D. y are the three Euler angles used to 


determine the relationship between (xo, y», zo) and (x, y, z)) 


获 规 则 采用 三 维 的 偏心 正八 面体 算法 或 二 维 
偏心 正方 算法 ,其 基本 原理 可 以 表述 为 : 若 某 一 优 
先生 长 方向 为 (a, B, 力 的 唱 粒 形 核 , 根据 其 大 计算 其 偏 
心 正 八 面体 顶点 到 其 中 心 的 距离 (0D, 其 计算 公式 为 : 
Lt + Ai) = L(t) + 42 Ax* Af. (18) 
其 中 , Ax 为 网 格 尺 寸 . 当 L(D) 求 得 以 后 , 就 可 以 知道 
偏心 正八 面体 的 6 个 顶点 在 局 部 坐标 系 Cco yo, zo) P 
的 坐标 位 置 , 然后 利用 式 (16) 就 可 以 计算 得 到 6 个 顶 
点 在 世界 坐标 系 (x, y, z) 中 的 坐标 位 置 . 若 在 某 一 时 
刻 正 八 面体 顶点 碰 触 到 周围 某 液 相 单元 时 , 该 液 相 
单元 就 会 被 捕获 成 为 界面 单元 , 并 赋予 优 先生 长 角 
度 值 (a, B, y). 在 后 续 的 生长 过 程 中 , 该 单元 也 会 依 
据 自 身 的 大 来 计算 它 的 正八 面体 顶点 坐标 位 置 , 从 
而 来 实现 对 周围 液 相 单元 的 捕获 . 
3 实验 研究 
3.1 实验 方法 
本 工作 采用 主要 成 分 为 AL-7Si-0.36Mg (质量 
数 , %) 的 铸造 铝 合金 , 熔炼 温度 为 720~740 C, 经 过 
Ar 气 除 气 处 理 , A1-5Ti-1B 细 化 处 理 以 及 Al-10Sr 变 质 
处 理 . 为 了 获得 不 同 的 凝固 条 件 , 将 处 理 后 的 金属 液 
浇注 到 具有 阶梯 形状 的 石英 砂 砂 型 型 腔 当中 , 同时 为 
了 记录 不 同 阶 梯 的 凝固 过 程 , 在 每 个 阶梯 的 中 心 位 置 
都 放置 了 一 个 K- 型 热电 偶 , 并 连接 到 温度 记录 仪 记录 
每 个 测 点 的 冷却 曲线 . 阶梯 铸件 的 形状 尺寸 以 及 热电 
偶 的 测 温 位 置 如 图 4 所 示 , 其 中 实心 点 处 就 是 热电 侦 
的 测 温 位 置 , 同时 把 热电 偶 的 测 温 位 置 按照 阶梯 的 厚 


$ 


度 从 厚 到 薄 依 次 从 No.1~No.7 进 行 标记 . 

为 了 分 析 不 同 凝固 条 件 下 的 组 织 形 貌 及 晶 粒 
尺寸 大 小 , 在 每 个 测 温 点 处 分 别 截取 直径 为 8 mm 
的 圆 形 截面 试 样 (No.7 处 的 试 样 直径 为 6 mm), 通过 
标准 的 制 样 工序 获得 金 相 试 样 . 本 工作 采用 阳极 履 
膜 和 偏振 光 金 相 技 术 来 区 分 不 同 的 晶 粒 . 阳极 履 膜 
剂 采 用 的 是 38 mL HSO,+43 mL HPO,+19 mL H,O 
配制 成 的 混合 溶液 , 电压 为 15 V, 电流 控制 在 0.05~ 
0.10 A 范围 内 , 履 膜 时 间 约 40~60 s. 78 RS Jes EI] ACIE 
在 Axio Scope.Al 型 且 带 有 和 补偿 器 的 金 相 显微镜 
(OM) 下 进行 观察 . 本 工作 采用 面积 法 来 测定 面 晶 粒 
密度 N, 其 公式 为 : 

Ni+0.5N,+1 


i MEET MECEEE (19) 


式 中 , NR ERA 内 完整 的 晶 粒 个 数 , 入 ;为 区 域 4 
边缘 上 的 晶 粒 数目 , 然后 根据 Poisson-Voronoi 模 型 站 
N, 20.568(N.)" , 计算 得 到 体 晶 粒 密 度 和 Nv. 二 次 枝 唱 
臂 间距 采用 截 线 法 来 测定 . 

3.2 实验 结果 与 分 析 

图 5 是 不 同 测 温 点 (No.1~No.7) 的 冷却 曲线 . 可 
以 看 到 , 该 合金 的 凝固 过 程 主要 包括 2 部 分 , 首先 是 
Q-Al 形 核 生 长 , 并 且 明 显 可 以 看 到 在 凝固 初期 , 由 于 
Q-Al 结 晶 潜 热 的 释放 , 在 温度 达到 Ti 之 后 , 会 出 现 
温度 的 回升 (再 辉 现象 ). 当 温度 达到 共 唱 温度 时 , 出 
现 (Al+SiD 共 晶 反 应 , 并 且 由 于 大 量 凝 固 潜 热 的 释放 ， 
温度 降低 速率 明显 减 慢 . 从 图 5 可 以 明显 看 到 , 随 着 阶 
梯 厚 度 从 No.1 的 40mm 减 小 至 No.7 的 4 mm, Tmin K 
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Thermocouple 


图 4 实验 用 阶梯 铸件 几何 尺寸 示意 图 


Fig.4 Schematic of dimensions of the step shape casting used to study the effects of cooling rate (unit: mm) 
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大 形 核 过 冷 度 AT,, H m JA BE T. 平均 冷却 速率 


700 

e li R, 晶 粒 密度 , 二 次 枝 唱 臂 间距 等 参数 . 从 表 1 中 可 以 
pr | 看 到 , 随 着 阶梯 厚度 变 小 , 熔 体 的 散热 作用 增强 , 晶 
2] . 粒 密度 从 增 大 , FRI Ye S AE. 
1 600 | f 为 了 求 得 式 (6) 当 中 的 参数 Nm 和 w 需要 利用 表 
P sl POAG d gd 1 中 的 AT, 和 ,数据 进行 曲线 拟 合 . 首先 对 式 (6) 两 

se0| 边 取 对 数 , 该 方程 可 以 简化 为 In(NW) 相 对 于 AT" 的 

9 ! 一 次 函数 形式 : 

0 100 200 300 400 500 
Timers ln(VJ=ln(V)-aA7 (20) 


SAl-7Si-0.36Mg 合 金 不 同 阶梯 处 测 得 的 冷却 曲线 
Fig.5 Cooling curves of different steps (No.1~No.7) with 


various thicknesses of Al-7Si-0.36Mg alloy sample 
casting ( Inset show the enlarged view of the nucle- 


ation interval on the cooling curves) 


次 降低 , 由 No.1 85 611.3 "C KEIR EI No.7 Rf) 595 "C, 
Za 的 降低 导致 形 核 区 间 增 大 , 从 而 使 得 熔 体 中 形成 
更 多 的 有 效 核 心 , 形 核 密度 增 大 . 

图 6 是 测 温 点 No.1~No.7 处 试 样 的 微观 组 织 照 
F. 由 于 三 维 空间 枝 晶 结构 的 复杂 性 , 二 维 平面 上 观 
察 到 的 枝 晶 并 不 是 很 有 规则 , 但 基本 可 以 体现 a- AT 
枝 晶 的 平面 四 重 对 称 结构 . 此 外 , 从 图 6 中 可 以 看 到 ， 
在 较 低 的 冷却 速率 条 件 下 , 由 于 形 核 数目 较 少 , 形 核 
位 置 不 均匀 , 使 得 同一 试 样 中 最 终 的 枝 唱 尺寸 差别 
较 大 , 表现 为 在 形 核 密度 高 的 位 置 枝 晶 尺寸 较 小 , 而 
在 形 核 密度 低 的 位 置 , 枝 晶 的 生长 不 受 限制 , 枝 晶 尺 
寸 较 大 . 随 着 冷却 速率 的 提高 , 形 核 密度 增 大 , 形 核 
核心 位 置 分 布 趋 于 均匀 , 最 终 的 枝 晶 尺 寸 差 别 不 大 . 
为 了 定量 表征 不 同 测 温 点 处 的 凝固 参数 和 组 织 特 


图 7 为 In(VW 与 A7 呈 的 函数 关系 . 从 图 中 可 以 
看 到 , No.1-No.6 If] In( Ny EDS T AT, 2 As E 2E TE 
系 , 但 是 No.7 的 数据 偏差 比较 大 . 其 主要 原因 是 测 
点 No.7 的 位 置 很 薄 (4 mm), 铸 型 壁面 引入 的 异 质 核 
心 对 形 核 过 程 影响 很 大 , 已 经 不 能 认为 TiB; 是 唯一 
的 形 核 核心 , 所 以 形 核 假设 条 件 (3) 不 再 满足 . 其 次 
是 该 处 的 散热 效果 很 强 , 难以 得 到 准确 的 测 温 数据 
(图 5), 可 能 会 导致 表 1 中 的 Ts 误差 较 大 . 所 以 , 本 
TE 3E No. 1-No.6 H In(NV) RA T, C ETT 
WE, 得 到 的 函数 关系 式 如 图 7 中 所 示 , 从 而 可 以 得 
到 Ny—e'75-2382.7 cm”, a=7.747 C, 代入 式 (6), 获得 
的 形 核 模型 为 : 


N, =2382.7 exp A) (21) 
依据 类 似 的 方法 , 也 可 以 得 到 面 形 核 密度 函数 : 
N, = 260.1 exp(- A (22) 


由 于 没有 考虑 No.7 数 据点 的 作用 , 本 工作 获得 
的 形 核 函 数 对 于 No.7 处 的 凝固 条 件 并 不 适用 , 在 后 
续 的 工作 中 还 需 进一步 完善 . 采用 该 模型 进行 形 核 


征 , 表 1 统 计 了 每 一 测 温 点 处 的 最 小 形 核 温度 Tuus, 最 


密度 预测 前 , 必须 先 要 通过 实测 或 是 模拟 得 到 铸件 
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Table 1 Experimental measured of solidification parameters and microstructure data for Al-7Si-0.36Mg alloy 
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6 测 温 点 No.1~No.7 的 微观 组 织 形 貌 和 晶 粒 尺寸 

Fig.6 OM images under polarized light showing the dendrite microstructure and grain size corresponding to the positions 
of No.1 (a), No.2 (b), No.3 (c), No.4 (d), No.5 (e), No.6 (f) and No.7 (g) (No.1-No.6: diameter 8 mm, No.7: diame- 


3 1 实验 测 得 的 凝固 参数 及 


多 核 密度 等 微观 组 织 数据 


No. 人 Atís R/(C:s)  N/cm?^ MNy/cm°  Aj/pum 
l 611.3 3.3 565.6 257.0 0.18 542 226.6 70.1 
2 611.0 3.6 563.1 208.0 0.23 64.1 291.5 66.8 
3 610.2 4.4 562.0 134.0 0.36 80.2 407.9 60.7 
4 608.8 5.8 560.8 83.0 0.58 101.6 581.6 53.2 
5 606.0 8.6 559.8 50.0 0.92 135.6 896.8 43.0 
6 604.1 10.6 558.7 30.0 1.51 172.1 1282.3 36.7 
7 595.0 19.6 557.0 7.5 5.00 418.2 4854.3 24.9 


Note: 7,,4—minimum nucleation temperature; AT,,—maximum nucleation undercooling temperature; 


T.—eutectic nucleation temperature; At—solidification time from Tmin to R; R.—cooling rate, 


calculated by R=(Tnn— T°)/At; Ns and Ny —surface and volume grain densities, respectively; 


secondary dendrite arm spacing. The liquidus temperature Ti, of this alloy is 614.6 'C, and 


is calculated by phase diagram software of Pandat 
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4 | 1 1 1 1 1 
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 
AT? [C? 
图 7Al-7Si-0.36Mg 合 金 的 InNv 与 AT,"' 的 函数 关系 
Fig.7 Linear relationship between grain density InNv and 


maximum nucleation undercooling AT; ' in the Al- 
781-0.36Mg alloy 
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9 某 一 点 处 相应 的 冷却 曲线 , 测 得 最 大 形 核 过 冷 度 
= AT m, 从 而 得 到 该 点 处 的 形 核 密 度 . 
"I INM =7.776-7.741A 7" 4 模拟 结果 与 分 析 

Tn . 4.4 二 维 枝 晶 生长 模拟 
0 本 工作 利用 建立 的 形 核 模型 和 三 元 枝 品 生长 
三 6f CA 算法 ,针对 No.1~No.7 测 温 点 处 的 微观 组 织 开展 
5. 相应 的 数值 模拟 . 为 了 体现 所 建立 的 形 核 模型 对 于 
该 合金 形 核 过 程 的 适用 性 , No.1~No.6 测 温 点 处 的 


形 核 数 目 由 式 (22) 计 算得 到 , No.7 处 根据 表 1 中 实 
测 值 代入 , 并 且 模 拟 所 需要 的 温度 值 通过 查询 实际 
测量 的 冷却 曲线 获得 . 模拟 所 采用 的 Al-7Si-0.36Mg 
合金 的 热 物 性 参数 如 表 2 中 所 示 . 模拟 所 采用 网 格 
数目 为 1800x1800, 网 格 尺寸 为 4 um, 梳 晶 的 优先 生 
长 角度 随机 赋予 .图 8 给 出 了 No.6 测 温 点 处 的 枝 晶 


表 2 Al-7Si-0.36Mg 合 金 模拟 所 需 的 参数 
Table 2 Parameters of Al-7Si-0.36Mg alloy used in the simulations 


Definition and symbol Value Unit 
Initial compositions, wi, wy, ws 77.0, wy, 70.36 96 
Liquidus temperature, T; (ws, wy) Calculated K 

Liquidus slope, mz, my, Calculated K-%" 

Partition coefficient, ks, kys Calculated 
Gibbs-Thomson coefficient, I 2.4x107 K.m 
Anisotropy coefficient, € 0.03 

Diffusion coefficient, D; Calculated m^.s' 

Timestep, 6 Ax'/6 max(D') s 


Note: “Calculated” means that the value is obtained through coupling with 


thermodynamic/kinetic/equilibrium phase diagram calculation databases 


Mass fraction of Si / % Mass fraction of Si / 96 


Mass fraction of Si / % 


Color online 
8 模拟 得 到 No.6 测 点 处 的 微观 组 织 演 化 和 Si 组 元 溶质 场 变化 过 程 
Fig.8 Evolution of simulated dendrite morphologies and the distributions of Si solute solidified at the position of No.6 for 
solidification time {=5.5 s (a), 9.7 s (b) and 20.1 s (c) 
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演化 过 程 和 溶质 场 分 布 . 从 图 中 可 以 看 到 , 在 凝固 
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凝固 组 织 . 可 以 看 到 , 随 着 冷却 速率 的 增 大 , 形 核 数 


初期 , 由 于 枝 唱 生长 较 慢 , 界面 前 沿 溶质 富 集 程度 


BUR, 枝 晶 以 四 重 对 称 的 形 


自由 生长 , 并 且 没 有 出 现 明显 的 枝 唱 臂 分 文 现象 (图 


貌 沿 着 各 自 的 晶体 取向 


8a). 随 着 温度 的 降低 , 枝 唱 生长 速率 v 增 大 (veA7”)， 


界面 前 沿 溶质 富 集 程度 加 
从 而 使 固 / 液 界面 形态 失 稳 


ul, 导致 成 分 过 冷 增 大 ， 
, 出 现 明 显 的 分 又 现象 ， 


二 次 枝 唱 臂 高 度 发 达 ( 图 8b). 此 外 , 随 着 枝 唱 不 断 长 
K, 相 邻 校 晶 间 的 生长 受到 相互 阻碍 . 这 种 阻碍 作 
用 也 主要 以 溶质 场 变 化 的 形式 体现 出 来 . 当 两 相 邻 


枝 品 相互 靠近 时 , 它们 之 间 


排出 的 溶质 由 于 没有 较 


好 的 扩散 通道 而 不 能 及 时 带 走 , 从 而 在 界面 前 沿 富 


集 形成 高 浓度 液 相 区 , 导致 


该 处 的 校 晶 生长 速率 减 


Te, 使 得 枝 晶 呈现 为 不 对 称 生长 (图 8b). 在 凝固 末 


期 , 剩余 液 相 区 的 溶质 高 度 富 集 , 枝 唱 尖端 的 生长 


受到 抑制 , 此 时 极品 的 进 一 


步 生 长 主要 通过 枝 品 辟 


的 生长 、 粗 化 ` 熔 合 等 方式 来 实现 , 而 二 次 校 晶 辟 间 
距 会 自动 的 调整 到 一 个 合适 的 范围 之 内 (图 8c). 
图 9 显示 的 是 不 同 凝固 条 件 下 模拟 得 到 的 最 终 


目 增多 , 晶 粒 尺寸 减 小 , 二 次 枝 晶 臂 不 断 细 化 , 模拟 
的 结果 与 图 6b, d, £, g 中 的 实验 结果 相符 合 . 为 了 定 
量 地 说 明 模 拟 结果 的 准确 性 , 图 10 给 出 了 不 同 测 温 
点 处 的 二 次 枝 晶 臂 间距 的 预测 值 和 实验 值 . 从 图 中 
可 以 看 到 , 模拟 结果 和 实验 结果 吻合 良好 , 两 者 的 误 
差 值 基本 控制 在 4% 以 内 . 在 文献 [7,30] 中 , 为 了 体现 
冷却 速率 对 枝 唱 臂 细 化 程度 的 影响 , 把 二 次 枝 晶 臂 
间距 与 平均 冷却 速率 及 通过 相应 的 索 函 数 关系 式 
=aRe* 关 联 起 来 , 其 中 a 和 4b 为 系数 . 本 工作 也 通过 
该 种 方式 对 模拟 得 到 的 数据 进行 了 拟 合 , 得 到 的 函 
TUS RUN A;742.8R. Y, 这 与 文献 [30] 中 通过 大 量 
实验 结果 拟 合 得 到 的 关系 式 =39.4R。"” 非常 接 
近 . 上 述 实 验 和 模拟 的 对 比 体现 了 该 模型 在 模拟 AL- 
7Si-0.36Mg 合金 枝 晶 生长 方面 的 准确 性 , 并 且 证 实 
该 模型 可 以 用 于 该 合金 二 次 枝 唱 臂 间距 的 预测 , 为 
后 续 的 性 能 预测 提供 了 前 提 条 件 . 
42 三 维 术 晶 生长 模拟 

由 于 枝 晶 的 实际 生长 过 程 是 三 维 空间 结构 , 二 


Colo 
9 模拟 获得 的 No.2, No.4, No. 
Fig.9 Simulated final dendrite microstructures of the positions No.2 (a), No.4 (b), No.6 (c) and No.7 (d) (Different colors 


represent different dendrites) 
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6 和 No.7 测 温 点 处 的 最 终 凝 固 组 织 
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维 平面 枝 唱 的 模拟 (图 多 很 难说 明 实际 观察 到 的 极 “， 演变 过 程 进行 模拟 是 非常 必要 的 , 目前 这 也 是 枝 晶 
其 复杂 的 枝 晶 形 貌 (图 6), 因此 从 三 维 角度 对 枝 晶 的 。 ”模拟 的 一 大 发 展 趋势 . 实际 的 凝固 过 程 涉及 到 多 枝 
晶 的 生长 , 并 且 相 邻 枝 晶 间 的 生长 会 受到 相互 的 阻 

Hw 碍 , 同时 还 会 受到 枝 晶 形 核 位 置 、 形 核 密度 、 枝 晶 优 
“| à Experimental —— 先生 长 取向 等 因素 的 影响 , 目前 已 有 学 者 "结合 X 
EE 射线 断层 扫描 技术 和 相 场 模拟 方法 开展 了 枝 晶 三 
维 结构 的 相关 研究 , 从 而 来 说 明 实际 观察 到 的 枝 唱 


80 


60- 


E 
x 50t 
» "I 14.242.8R 三 维 结构 的 复杂 性 和 多 样 性 . 本 工作 利用 建立 的 
CA 模型 对 Al-7Si-0.36Mg 合 金 的 三 维 枝 晶 生长 过 程 
Me 和 枝 晶 形 貌 开展 数值 模拟 . 模拟 采用 的 网 格 量 为 
: ] 400x400x400, 网 格 尺寸 为 5 pm (计算 域 大 小 2 mmx 
R.I CCS’) 2 mmx2 mm). 本 工作 针对 No.7 测 温 位 置 开展 三 维 


枝 品 生长 模拟 研究 , 其 中 计算 域 中 的 形 核 数 目 根据 
表 1 中 实测 的 形 核 密度 数据 得 到 , 且 每 个 晶 核 都 随 
机 赋予 生长 取向 , 温度 数据 通过 读 取 实 测 的 温度 | 


10 不 同 测 温 点 位 置 的 三 次 枝 蝇 臂 间距 模拟 结果 和 实 
验 结果 对 比 
Fig.10 Comparison between the predicted secondary den- 


drite arm spacing and the experimental ones solidi- 线 来 获取 . 图 11 所 7R 的 是 No.7 位 置 凝 [zi] 条 件 下 的 三 
fied in different cooling conditions 维 多 TX HH 2H 织 演 变 过 程 . 从 图 中 可 以 看 到 ， 由 F 在 


Color online 
11 Al-7Si-0.36Mg 合 金 在 No.7 位 置 处 的 三 维 多 枝 唱 组 织 演变 过 程 


Fig.11 Simulated evolution of multi-equiaxed dendrites for Al-7Si-0.36Mg alloy solidified at the position of No.7 (f.—solid fraction) 
(3)—2.396,:-1.4s — (b)f-14.896,-24s — (c)f-25.199,1-33 s. (d)fi-54.896, 177.6 s 


Color online 
12 凝固 结束 时 二 维 截面 枝 晶 形 貌 和 实验 结果 对 比 


Fig.12 Simulated final dendrite morphologies on arbitrarily two dimensional sections (a) and the obtained experimental 


metallographic microstructure solidified at the position of No.7 (b) 
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凝固 初期 枝 晶 尺寸 较 小 , 相互 影响 较 弱 , 枝 晶 可 以 
在 过 冷 炊 体 中 沿 优先 生长 方向 自由 生长 (图 11a). 随 
着 温度 的 降低 , 过 冷 度 增 大 , 枝 晶 生长 速度 加 快 , 界 
面 稳定 性 减弱 , 一 次 梳 晶 主干 出 现 了 二 次 分 支 , 了 
且 由 于 受到 相 邻 枝 唱 或 是 计算 域 边界 的 影响 , FE 
开始 呈现 不 对 称 生 长 (图 11b 和 c). 在 凝固 末期 , 枝 
占 满 整个 计算 域 , 枝 晶 辟 不 断 粗 化 、 熔 合 , 使 得 最 
呈现 为 结构 复杂 、 形 貌 多 样 的 枝 晶 组 织 ( 图 11d). 
12 呈现 了 初生 相 凝 固 结束 时 二 维 截 面 枝 晶 形 貌 
和 实验 结果 比 对 . 从 图 中 可 以 看 到 , 二 次 枝 唱 臂 呈 
断 续 分 布 , 相对 于 二 维 模 拟 结果 (图 9), 三 维 条 件 下 
模拟 得 到 的 枝 晶 形 貌 与 实验 结果 更 加 有 具有 相似 性 . 
此 外 , 利用 截 线 法 对 图 12a 和 b 中 实验 和 模拟 得 到 的 
No.7 位 置 处 的 二 次 枝 唱 臂 间距 进行 统计 , 2 者 分 别 
为 28.9 和 24.9 um. 上 述 的 模拟 结果 充分 体现 了 本 
模型 在 模拟 三 维 枝 晶 生长 方面 的 可 行 性 . 
5 结论 

(1) 通过 分 析 阶 梯 件 不 同位 置 处 的 冷却 曲线 , 3 
基于 合理 的 形 核 条 件 假 设 , 建立 了 适用 于 普通 凝固 条 
件 下 Al-7Si-Mg 合金 等 轴 唱 凝固 的 形 核 模型 . 该 模型 
中 形 核 密度 Ns 随 最 大 过 冷 度 的 倒数 AT' 呈 指数 性 变 
化 , 并 且 该 模型 中 的 相关 参数 都 有 明确 的 物理 意义 . 

(2) 通过 与 热力 学 、 动 力学 、 平 衡 相 图 数据 库 相 
BE, 并 利用 空间 坐标 变化 等 算法 , 建立 了 适用 于 
三 元 铝 合金 二 维 、 三 维 枝 晶 生长 的 CA TUS. 利用 建 
立 的 形 核 和 生长 模型 , 针对 Al-7Si-0.36Mg 合金 开展 
了 二 维 枝 晶 生长 过 程 数 值 模拟 , 分 析 了 不 同 凝固 条 
牛 下 的 枝 晶 演化 和 形 貌 特征 , 并 且 模拟 得 到 的 二 次 
枝 晶 臂 间距 与 实验 吻合 良好 , 体现 了 该 模型 在 模拟 
三 元 铝 合金 枝 唱 生 长 方面 的 适用 性 . 

(3) 模拟 了 Al-7Si-0.36Mg 合金 三 维 多 枝 唱 生 长 
过 程 , 并 在 二 维和 截面 上 与 实验 结果 进行 对 比 , 初步 
分 析 了 枝 唱 三 维 空间 结构 复杂 性 和 多 样 性 , 体现 了 
本 模型 在 模拟 三 维 枝 晶 生长 方面 的 可 行 性 . 
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